Orbital state and magnetic properties of LiV204 
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LiV2 04 is one of the most puzzling compounds among transition metal oxides because of its heavy 
fermion like behavior at low temperatures. In this paper we present results for the orbital state and 
magnetic properties of LiV204 obtained from a combination of density functional theory within the 
local density approximation and dynamical mean- field theory (DMFT). The DMFT equations are 
solved by quantum Monte Carlo simulations. The trigonal crystal field splits the V orbit als such 
that the ai^ and e^ orbitals cross the Fermi level, with the former being slightly lower in energy and 
narrower in bandwidth. In this situation, the d-d Coulomb interaction leads to an almost localization 
of one electron per V ion in the ai^ orbital, while the e^ orbitals form relatively broad bands with 
1/8 filling. The theoretical high-temperature paramagnetic susceptibility x{T) follows a Curie- Weiss 
law with an effective paramagnetic moment pe//=l-65 in agreement with the experimental results. 



I. INTRODUCTION 

Heavy-fermion (HF) materials are typically intermetal- 
lic compounds containing Ce, U, or Yb atoms. They 
are characterized by extraordinarily strongly|-.mnormal- 
ized effective masses, m* ^ 100 - 1000 rngP'oB as in- 
ferred from the electronic specific heat coefficient 7(T) = 
Ce{T)/T at low temperature. They also show an appar- 
ent local moment paramagnetic behavior with a strongly 
enhanced spin susceptibility x low temperatures. The 
discovery by Kondo et alu of the HF behavior in LiV2 04 
with a Kondo or spin fluctuation temperature T^ ~ 
28 K has significant importance, because this is the 
first c?-electron system that shows HF characteristics. 
Kondo et al. reported a large electronic specific heat 
coefficient 7 ~ 0.42J/(mol K^) at 1 K, which is much 
larger than those of other metallic transition metalj-com- 
pounds like, e.g., Yi_^Sc^Mn2 k < 0.2 J/(mol K^)p and 
[7 <0.07 J/(mol K^)].ll Also a crossover from 
the local moment to a renormalized Fermi-liquid behav- 
ior was obsewed with decreasing temperature. Recently 
Urano et aln reported that the electrical resistivity p 
of single crystals exhibits a T^ temperature dependence 
P = pQ-^ with an enormous A_which like in conven- 
tional HF systems scales with 7^.@ In the temperature 
range 50-1000 K the experimental magnetic susceptibil- 
ity follows a Curie- Weiss law with a negativ Curie- Weiss 
temperature, which indicates a weak antiferromagnetic 
(AF) V-V spin interaction (see later in the text). No 



magnetic ordering was observed down to 0.02 K.I 

These unexpected phenomena entailed numerous ex- 
periments, which confirmed the HF behavior of Liy204 
in a variety of physical quantities: Johnston et aln car- 
ried out specific heat and thermal expansion measure- 



ments. Kondo et a/.0 described the synthesis, char- 
acterization and magpetic susceptibility versus temper- 
ature. Onoda et alx3 explored spin fluctuations and 
transport in the spinel systems LiajMgi_a^V2 04 and 
Lia;Zni_a;V2 04 through measurements of x-ray diffrac- 
tion, electrical resistivity, thermoelectric power, magne- 
tization, and nuclear magnetic resonance (NMR). The 
electron-spin resonance (ESR) and the magnetic sus- 
ceptibility in purcp^uid doped LiV204 were measured 
by Lohmann et al.xB Photoemission studies of the hole 
doped Mott iasulator Lii_a^Zna^V2 04 were done by Fu- 
jimori et al.M A series of ^Li-=]NMR experiments were 
carried out by Fujiwara et al.t3 for LiV2 04 and also 
for Lii_a;Zna^V204.LJ The Knight shift, spin suscepti- 
bility and relaxation times were determined from ^Li^ 
NMR experiments by Mahajan et al.Lj Krimmel et al$3 
presented results for the magnetic relaxation of LiV204 
obtained by meeps of quasi-elastic neutron scatter- 
ing. Trinkl et al&B investigated spin-glass behavior in 
Lii_a;Zna;V2 04. Urano et aln experimentally observed 
results for C(T), x(T), resistivity p{T) and Hall co- 
efficient Erp(T) for single crystal samples. Recently, 
Lee et a/.E2l performed inelastic neutron scattering mea- 
surements and Fujiwara et alE3 studied the spin dynam- 
ics under high pressure. A review of various experi- 
ments and tbeoretical research is collected in the work 
of Johnston.EHl 

Let us summarize the experimental efforts concern- 
ing the measurements of the spin susceptibility x(T). 
As mentioned before, a temperature dependent spin sus- 
ceptibility x(T) is observed in the temperature range 
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x(T)=xo+C/(T-6>). However, depending on the quality 
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of samples and experimental techniques, the values of the 
Curie constant C obtained from fitting the experimeiit|ai 
data are in the range from 0.329-0.468 cm^K/mol.BEl 
More important and of direct physical importance is the 
value of the effective paramagnetic moment Pef / , which is 
defined via the Curie constant asp^^j := 3CkB/{NA/J^%)' 
The extmjjted value Pe// lies in the range between 1.62 
and 2.1.E3Ej The Curie- Weiss temperature obtained in 
those experiments is also slightly different and encom- 
passes the interval = —20 ... — 60 K, which indi- 
cates a weak antiferromagnetic (AF) V-V spin interac- 
tion. If one further assumes that the effective spin value 
is S=l/2, an experimental estimation of the Lande fac- 
tor g'^ =j)lff/[S{S-\- 1)] leads to values between 1.87and 
2.23.00 

From a theoretical point of view, LiV204 has been in- 
tensively studied by standard band structure calculations 
by various groups usi]m,different implementations of den- 
sity functional theoryuiL within the local density approx- 
imation (DFT/LDA).I23 Anisimov et alx^ investigated 
the possibility of localization of the d-states in LiV2Q4| 
within the linearized muffin-tin orbital^XLMTO) basisBll 
supplemented by the LDA+U method.LJ The elGdjronic 
structure was furthermore studied by Eyert et alSB with 
the scalar-relativistic augmented spherical wave (ASW) 
basis,LZl while Matsuno et a/J23 jused a full-potential, 
scalar-relativistic implementationEj of j-the linear aug- 
mented plane wave (LAPW) approached for the band 
structure calculation and fp;thermore applied a simple 
tight-binding (TB) mockirti^ingh et alJt3 also used 
the full potential LAPWBc^S foij-ealculating the band 
structure and a TB-LMTO methodo to analyze the band 
symmetry. 

A relative comparison of LiV204 and lithium titanai^ 
(LiTi2 04) with original HF systems was done by Varmao 
to provide some qualitative understanding of these com- 
pounds. The two-band Hubbard model in the slave- 
boson mean-field approximation was applied to LiV204 
by Kusunose et alxll to investigate the evolution of bands 
due to the Coulomb interactions— Hopkinson et al. pre- 
sented a simple two-band modeled and a t-Xjpiodel with 
a strong Hund's coupling for the (i-electronst£l to find ev- 
idence for a two- stage screening in LiV2 04. Fujimoto in- 
vestigated the Hubbard chains network model on corner- 
sharing tetrahedra as a po^ible microscopic model for 
the HF behavior in LiV204.L3 The competition between 
the Kondo effect and frustrating exchange interactions 
in a Kondo-lattice model within a large- A/" approach for 
the spin liquid together with dynamical mean-field the- 
ory were studied by Burdin et a/. .Ell 

Despite all this theoretical effort, there does not yet 
exist an undisputed microscopic explanation of the HF 
behavior of LiV204 at low temperatures. It has been 
attempted to explain the low-temperature properties of 
LiV204 by a mechanism analogous to one for systems 
with local moments on a pyrochlor lattice, which are frus- 
trated with respect to the antiferromagnetic interactions. 
A novel attempt to provide a microscopic model based on 



material properties of LiV2 04, which does not evoke the 
idea of frustration was suggested by Anisimov et al.x^ 
Their basic idea is a separation of the electrons on par- 
tially filled t2^ orbitals into localized ones, forming lo- 
cal moments, and delocalized ones, producing a partially 
filled metallic band. The hybridization between those 
two subsets of electrons, like in /-electrons materials, can 
give rise to heavy- fermion effects. However, in this wofk 
the conclusions were based on an LDA+U calculation,^^ 
which is essentially a static mean-field approximation and 
surely too crude to give a proper description of LiV204. 

The present work reports on the investigation of the 
orbital state and the magnetic properties of LiV204 us- 
ing the novel LDA+DMFT(QMC) approach. This mate- 
rial specific theory extracts informations about the non- 
interacting band structure from a DFT/LDA calculation, 
while treating the local Coulomb interactions via the 
dynamical mean- field theory (DMFT), which is a well- 
established non-perturbative approach to study localiza- 
tion effects in strongly correlated materials. The result- 
ing DMFT equations are then solved by quantum Monte- 
Carlo simulations. The main goal of this study is to fur- 
ther explore the idea of a separation of the d-electrons in 
LiV204 into two subsets, localized and itinerant as pro- 
posed in Ref. The rest of the paper is divided into two 
parts: In the first we discuss standard DFT/LDA results 
(section II) and then present our LDA+DMFT(QMC) 
results in the second (section p^ ). A short summary 
concludes the paper. 



II. LIV2O4: DFT/LDA RESULTS 
A. Crystal structure and d orbital splitting 

LiV204 has the fee normal- spinel structure with non- 
symmorphic space group Fd3m aiid was first synthesized 
by Renter and Jaskowsky in 1960.E3 The Li ions are tetra- 
hedrally coordinated by oxygens, while the V sites are 
surrounded by a slightly distorted edge-shared octahedral 
array of oxygens (Fig. |l|). The corresponding unit cell of 
the face-centered cubic lattice contains two LiV204 for- 
mula units (14 atoms) with four V atoms, which form a 
tetrahedron (Fig. ||). The LiV2 substructure is the same 
as the C15 structure AB2, where-the local moments at 
the B sites are highly frustrated.ES The observed lattice 
constant of LiV2 04 is 8.22672 A at 4 kM The eight 
oxygen atoms in the primitive cell are situated at the 
32 e- type sites, at positions which are determined by the 
internal-position parameter x=0.2611 in units of the lat- 
tice constant. 

The total space group of the crystal is cubic but the 
local point group symmetry of the V ion crystallographic 
position is trigonal Dsd- The different trigonal axes of ev- 
ery V atom in the unit cell are directed towards the center 
of the tetrahedron (Fig. |^). Since the formal oxidation 
state of V is non-integer V^"^^ (configuration d^'^) and 
the V ions are crystallographically equivalent, LiV204 



3 





FIG. 1: The normal- spinel crystal structure of LiV2 04 is 
formed by the oxygen edge-shared octahedra with V atoms 
at the centers. V: large dark ions; O: small dark ions. 
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FIG. 3: Schematic picture of the electronic level distribu- 
tion caused by the strong p-d hybridization in the trigonally 
distorted octahedral crystal field. The notations i2g and e^ 
correspond to the irreducible representations of the cubic Oh 
group, aip and e^ to the irreducible representations of the 
Dsd group. 

(Fig. H). However, this splitting is not Igjge enough to 
separate the t2g band into two subbands.H 



B. DFT/LDA band structure 
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FIG. 2: Tetrahedron formed by four V atoms (large dark 
spheres) in the spinel unit cell and the trigonal axes for each 
V atom. 



Based on DFT/LDA0'i within the LMTO method,i 

we performed first-principle calculations of the electronic 
structure of LiV2 04. The radii of the muffin-tin spheres 
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must be metallic as it is observed.ES 

The octahedral crystal field at the V sites in the spinel 
structure splits the V 3d bands into three degenerate and 
partially filled t2g bands and two empty e^ bands. The 
Fermi level lies within the t2g complex, thus the transport 
properties of LiV204 are solely associated with the t2g 
bands. The trigonal symmetry of the V ions splits the 
cubic t2g levels into one ai^ and two degenerate e^ levels 



FIG. 4: DOS of LiV204 calculated with the LDA-LMTO 
method. Upper figure: Y-Sd (full line) and 0-2p (dashed 
line) DOS; lower figure: partial Y-3d{t2g) (full line) and V- 
3d{eg) (dashed line) DOS. The Fermi level corresponds to 
zero energy. 

were RLi=2.00 a.u., Rv=2.05 a.u., and Ro=1.67 a.u. 
The resulting densities of states (DOS) are shown in 
Figs. ^ and ||. In Fig. ^, there are three well-separated 
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sets of bands: completely filled 0-2p-bands, partially 
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FIG. 5: Partial Y-3d{t2g) DOS of LiV204 calculated with 
the LDA-LMTO method: aig (full hne) and e^ (dashed 
line) projected DOS. The center of gravity of the aig or- 
bital is Caig =0.14 eV; that of the e^ orbital is Ce^ =0.40 eV. 
The LD A- calculated occupancies of the aig and e^ orbitals, 
naig=0.60 and ne^=0.45, are nearly the same. The Fermi 
level corresponds to zero energy. 

filled t2g bands, and empty e^ bands. The bands in the 
energy range from -8 eV to -3 eV originate mainly from 
0-2p states and have only small admixtures from V-Sd 
states. The upper two groups of bands, which extend 
from -1.0 eV to 0.8 eV and from 2.3 eV to 3.2 eV, are 
predominantly derived from the V-3d states. Although 
additional small 0-2p contributions are apparent in this 
energy range, the p-d hybridization is much reduiped com- 
pared to other early transition metal oxides. E3 Due to 
the crystal field of the slightly distorted octahedral co- 
ordination of the V atoms by the oxygen atoms, a clear 
separation of the 3d i2g and e^ groups of bands is visible 
in Fig. ^. Whereas the former states appear exclusively 
around the Fermi energy, the e^ states prevail at higher 
energies. Contributions of the \ 3d states to the oxygen- 
derived bands originate almost exclusively from the e^ 
states, which form a-bonds and experience a strong over- 
lap with the 0-2p states. In contrast, the i2g orbitals, 
which give rise to tt bonds, yield only a negligible con- 
tribution in this energy range. In addition to p-d bonds, 
the t2g states experience strong cr-type overlap with the 
t2g orbitals at neighboring vanadium sites of the fee sub- 
lattice. Hence, these (i-states take part in two different 
types of bonding, namely a-type V-V and 7r-type V-0 
bonding, which leads to two different bandwidths for i2g 
and e^ states. Since both the metal-metal and metal- 
oxygen bonds are mediated through the same orbital, a 
simple analysis of the partial DOS would not allow to 
distinguish tlie different roles played by the i2g orbitals. 
Eyert et. alx^ used a local coordinate system with the Z- 
axis along the trigonal direction (111) pointing towards 
the center of the tetrahedron formed by V ions (see Fig. 
^) and plotted the V d^^z'^-r'^-, dxz +dyz and dx2_y2 -\- dxy 
partial DOS. While d^z'^-r'^ orbitals are of pure t2g char- 



acter, the other four orbitals comprise a mixture of t2g 
and eg states. We have calculated the partial DOS for 
diig and e^ orbitals, using the irreducible repr£isentations 
of the Dsd group according to Terakura et al.l^ with the 
following linear combinations of the t2g cubic harmonics: 
the diig orbital is given by {xy -\- xz -\- yz)/\^ and the 
two eg orbitals by {zx — yz) / \f2 and (%jz zx — 2xy) / \/6- 
These three particular linear combinations are valid if 
the coordinate axes X, Y and Z are directed along V- 
O bonds. If the Z direction is chosen along one of the 
trigonal axes described above, then the ai^ orbital is the 
3^2 _ ^2 orbital in the local coordinate system. 

The projected partial LDA DOS of the ai^ and the e^ 
orbitals are shown in Fig. ^. The bandwidth of the ai^ or- 
bital Waig = 1.35 eV is almost a factor of two smaller than 
the e^ bandwidth We^=2.05 eV. Nonetheless, we found 
within LDA for all three t2g derived orbitals nearly the 
same occupancies: naig=0.60 and ne7r=0.45. Almost all 
spectral weight of the ai^ orbital is concentrated around 
the Fermi level in the region from -0.2 eV to 0.3 eV. In 
contrast to the ai^ orbital the e^ DOS is flat at the Fermi 
level. The largest part of the spectral weight of the e^ 
orbitals is situated in the interval from 0.3 eV to 0.85 eV. 
Despite such a different spectral weight the ai^ and e^ 
bands are not completely separated in energy. However, 
there is a significant difference in the centers of gravity 
calculated from the corresponding partial DOS, which 
can be interpreted as a measure for the trigonal splitting 
of the t2g states. 

The trigonal splitting is much smaller than the band- 
width but has a great importance for the understanding 
of the physics of the LiV204 system in the presence of 
strong Coulomb interaction. The value and the sign of 
the trigonal splitting will determine the orbital in which 
the V 3d electrons should be localized when a strong 
Coulomb interaction, jaihich is larger than the bandwidth, 
is taken into account .l3 We found that this trigonal split- 
ting value is very sensitive to the accuracy of the band 
structure calculations. In order to increase the accuracy, 
the overlap between atomic spheres in our LMTO calcu- 
lation was set to zero and more empty spheres were intro- 
duced. For the following we define the trigonal splitting 
of the t2g states as the difference of the centers of gravity 
of the diig and the ej projected DOS. We find that the 
diig center of gravity is 0.26 eV lower than that of the e^ 
band (caig=0.14 eV, Ce^=0.40 eV; see Fig. H). One can 
thus conclude from our LDA calculations that the ai^ 
orbital is more favorable for the localization of electrons. 

The effective mass m* is known to be a measure for 
the Coulomb correlations and can be obtained from the 
electronic specific heat coefficient 7 at low temperatures. 
From the LDA-calculated electronic specific heat coef- 
ficient 7^^^, one can infer the ratio of the effective 
mass to the band mass m*/m5=7/7^^^. Here, the 
former is related to the LDA DOS at the Fermi level 
via 7^^^=7r^/c^A/'y^D(£^i?)/3, and 7 is the experimen- 
tal value of the electronic specific heat coefficient taken 
from.a For LiV204 we found m*/m5~ 25.8, which is in 
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good agreement with previous results. ESI Such a huge en- 
hancement of the quasi particle mass m* is a strong 
evidence that Coulomb correlations are important in 
LiV204 and have to be taken into account in order to 
describe the physics of this system. 

III. A MICROSCOPIC THEORY FOR LIV2O4 
A. LDA+DMFT(QMC) scheme 

Based on the d = . ]lm[ the dynamical mean- 
field theory (DMFT)HBE3 was developed as a non- 
perturbative approach to describe strongly correlated 
electron systems. It permits the calculation of electronic 
spectra for systems with local electronic Coulomb cor- 
relations. The LDA+DMFT approach is ^rmeigei of 
the DFT/LDA and the DMFT techniques. BeI It com- 
bines the strength of the DFT/LDA, viz describing the 
weakly correlated part of the ab initio Hamiltonian, 
i.e., electrons in s- and p-orbitals and the long-range 
part of the Coulomb interaction, with the ability of the 
DMFT to treat electronic correlations induced by the 
local Coulomb interaction. In this paper we will only 
briefly discuss the relevant parts of the LDA+DMFT ap- 
proach and refer the reader to a recent report by Held et 
alx3 for more details. 

For a given material, one can extract a LDA Hamilto- 
nian ^^LDA supplement it the local Coulomb inter- 
actions: 

^ = ^LDA + U'^'^him^himi (1) 

m i 

Here, the index i enumerates the V sites, m denotes 
the individual t2^ orbitals, and a the spin. i^L^A is a 
one-particle Hamiltonian generated from the LDA band 
structure with an averaged Coulomb interaction sub- 
tracted to ^avoid a double counting of the Coulomnb 
inter act ion. EJ IJ is the local intra-orbital Coulomb re- 
pulsion and J the exchange interaction. The local inter- 
orbital Coulomb repulsion U' is then fixed by rotational 
invariance: U' = U — 2J . The actual values for U and 
U' can be obtained from an averaged Coulomb parame- 
ter U and Hund's exchange J, which can be calculated 
with LDA. The quantity U is related to the Coulomb 
parameters U and U' via 

j-j _ U + (iVorb - l)U' + (iVorb - IW - J) 

2iVorb-l ' 

where Norb is the number of interacting orbitals (Norb=3 
in our case). Since U and U' are not independent, the 
two values Cf-fwi J are sufficient to determine U from 
this relation.H'El 

The DMFT maps the lattice problem (|l|) onto an 
effective, self-consistent impurity problem. A reliable 



method to solve this (multi-band) quantum impurity 
problem is provided by quantum Monte Carlo (QMC) 
simulations,!^ which are combined with the maximum 
entropy methodtJ for the calculation of spectral functions 
from the imaginary time QMC data. This technique has 
been applied pta-ejaic properties of several transition 
metal oxides.ElOEl 

A computationally important simplification is due to 
the fact that in cubic spinel the i2g states do not mix 
with the e^ states. In this particular case the self-energy 
matrix T^rnaiz) is diagonal with respect to the orbital and 
spin indices. Under this condition the Green functions 
Gma{z) of the lattice problem can be expressed in the 
DMFT as Hilbert transform of the non-interacting DOS 

f (e) 
Gmcr{z) = / de . (3) 

This procedure avoids the rather cumbersome and prob- 
lematic /c-integration over the Brillouin zone by the an- 
alytical tetrahedron method.0 The LDA DOS N^\e) 
for different orbitals (different m) of the i2g states in 
LiV204 is displayed in Fig. ||. 

A particularly interesting quantity for LiV204 is the 
magnetic susceptibility x(^)- order to calculate x(^) 
we use the definition x{T) = lim^^o ^(^)/^, where 
M{T) is the magnetization due to the applied magnetic 
field B. Obviously, M{T) can be obtained directly from 
the QMC data via 

A^orb 

M(T)=AiB^((n™T)-(n„^i)) (4) 

m=l 

and 

{Uma) = G^„{t = 0+) . (5) 

Since QMC results are subject to both statistical and sys- 
tematic errors, a direct evaluation of M(T)/B for small 
B will inevitably lead to unpredictable scattering of the 
results for x{T). To avoid these problems, we perform 
calculations for M(T^ B) for a series of small fields B 
and extract x(^) from a least-squares fit as the slope of 
M(T, B) diS B ^ 0. 



B. Single-particle properties 

The DMFT calculations are based on the DFT/LDA- 
DOS N^{e) for the ai^ and e^ orbitals presented in Fig. ||. 
The total number of electrons in these three orbitals was 
fixed to n=1.5, in accordance with the +3.5 valency of 
V ions in LiV204. The d-d Coulomb interaction param- 
eter U=3.0 eV and exchange Coulomb interaction pa- 
rameter J=0.^-|eV were calculated^ by the constrained 
LDA method.EJ The temperature used in our QMC sim- 
ulations was approximately 750 K. While the scheme in 
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FIG. 6: Partial ai^ and DOS for LiV204 calculated with 
LDA+DMFT(QMC) for IJ = 3.0 eV and J = 0.8 eV. The 
non- interacting DOS N^(e) used in the LDA+DMFT calcula- 
tions are presented in Fig. g The LDA+DMFT(QMC) occu- 
pancies of the aip and e^ orbitals are Ua^^ =0.88 and =0.31. 
The Fermi level corresponds to zero energy. 



principle poses no restrictions on the temperature val- 
ues, the QMC code used presently limits our calculations 
to these rather higk temperatures because of computing 
power limitations]^ 

The partial DOS of ai^ and e^ orbitals obtained from 
the analytical continuation of the QMC results are shown 
in Fig. ^. In comparison with the non-interacting case 
(see Fig. a considerable transfer of spectral weight 
has taken place, especially the ai^ orbital appears to be 
strongly renormalized by Coulomb interactions. Further- 
more, in contrast to the LDA results, where the occupa- 
tion of the ai^ and ej orbitals is roughly the same, the 
DMFT yields n^,^ 0.9 and ^ 0.3. These numbers 
clearly demonstrate that one electron is nearly localized 
in the ai^ orbital, while the e^ states remain weakly cor- 
related and metallic with a filling close to 1/8. Nev- 
ertheless, the e^ orbitals are remarkably broadened by 
Coulomb interactions in comparison with the LDA pic- 
ture (fig. 1^). Slight differences between the two degener- 
ate e^ orbitals originate from the statistical nature of the 
QMC method and the following analytical continuation 
withh MEM. 

To discuss the origin of the structures in the DOS 
of Fig. |[ it is helpful to look at the spectrum of the 
atomic Hamiltonian and the corresponding positions of 
the one-particle excitations. The states corresponding 
to atomic occupancies up to two electrons are listed in 
Tab. I together with their energies. Since the Hubbard U 
is large and we have to accommodate 1.5 electrons, with 
one electron in the ai^ orbital, the ground state has to 
be a suitable mixture of two states, which both have a 
singly occupied ai^ orbital, but differ in the occupancy 
of the e^ orbitals. The inspection of all possibilities in 
Tab. I leaves only |cr)|0)|0), |cr)|cr)|0) and |cr)|0)|cr) as can- 
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FIG. 7: LDA+DMFT(QMC) calculated DOS for LiV204 
(solid line; [/==3.0 eV and J=0.8 eV) in comparison 
with non- interacting LDA ai^ and e^ DOS (dashed line). 
Please note that for a more convenient comparison the 
LDA+DMFT(QMC) DOS were magnified in intensity by a 
factor of 5. The Fermi level corresponds to zero energy. 
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TABLE I: Eigenstates and -energies of the atomic Hamilto- 
nian for total occupations less or equal two. ta^g denotes the 
one particle energy of the ai^ orbital and Ae ~ 0.26 eV the 
trigonal splitting between ai^ and e^ orbitals. 



didates. In order for these three to be (nearly) degener- 
ate, e^i^ + Ae + - J = or e^^^ = -V + J - Ae. 
Inserting the numbers for U' and J and the LDA- value 
for Ae leads to e^i^ ~ —3.1 eV. The possible one-particle 
excitations with respect to this ground state configura- 
tion and their energies can now be constructed easily, 
leading to the results in Tab. ||. It is important to note 
that the single-particle excitations of the ai^ orbital with 
cc^o < have two distinct contributions, namely one of the 
usual type "singly occupied unoccupied" , and a second 
that actually involves a doubly occupied state build of a 
mixture of ai^ and e^ states, i.e., a mixed- valent state. 
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TABLE II: Single-particle excitations for the atomic model 
filled with 1.5 electrons. For spin-degenerate processes only 
one representative is listed. 



The energy of this latter excitation is given directly by 
the trigonal splitting, i.e., this feature can also serve as 
a means to extract this number from photoemission ex- 
periments. 

That the excitations listed in Tab. || directly map to 
the peaks in the DOS Fig. || can be shown by using an- 
other technique to solve the DMFT equations, namely 
resolvent perturbation theory and NCA (see e.g. Ref. 68). 
This approach allows for a direct identification of differ- 
ent excitation channels (i.e., different initial particle num- 
bers) and a distinction between ai^ and e^ states. The 
result is shown in Fig. ^, where the black lines denote ai^ 
and the gray ones e^ single-particle excitations. Full lines 
stand for a singly occupied initial state, dashed lines for 
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FIG. 8: One-particle DOS separated into different contribu- 
tions. Black (thick) lines denote aig (e^) single-particle exci- 
tations. Full (dashed) lines stand for a singly occupied (dou- 
bly occupied) initial state, dot-dashed lines for a triply occu- 
pied initial state. The dotted vertical lines show the positions 
of the excitations listed in Tab. 

a doubly occupied initial state and dot-dashed lines for 
triply occupied initial state. Each of the dominant peaks 
correspond to one of the transitions listed in Tab. ||. 



For clarity, the different excitation energies are marked 
by dotted vertical lines. Apparently, all peaks in Fig. ^ 
have their counterparts in Fig. |[ explaining them as sim- 
ply due to the atomic multiplet structure. The shifts are 
typical renormalizations occurring in quantum impurity 
models. Especially the sharp peak at ~ —0.25 eV in 
the ai^-DOS is part of the lower Hubbard band, i.e., the 
Fermi energy is located at the upper edge of the lower 
Hubbard band, leading to the observed filling close to 1. 
These explanations show that no Kondo-resonances are 
seen in the LDA+DMFT(QMC) spectra at these elevated 
temperatures, consistent with experiment. Note also that 
the structures in the e^-band around uj ^ —1.5 eV are 
related to corresponding features in the LDA e^-band lo- 
cated roughly 1.5 eV below the main structure at 0.75 eV 
(see Fig. ^). The double peak structure at cc; > 0, on the 
other hand, consists of two different processes with exci- 
tation energies cjq = and coq = J, respectively. Finally, 
the peak around a; = 4 eV in the e^ DOS can be identi- 
fied as an excitation into a triply occupied state \a)\a)a) 
with energy U — 3 J. 

The orbital occupation, which confirnis-|the earlier con- 
jecture based on the LDA-hU approach,L3 can be readily 
understood from the previous analysis of the LDA band 
structure. There we found that the center of mass of 
the diig orbital is 0.26 eV lower than the corresponding 
value for the e^ orbitals. In the absence of Coulomb cor- 
relations the bandwidth is significantly larger than this 
energy difference, i.e., this splitting does not have any 
significant effect. However, with a Coulomb interaction 
U=3.0 eV, which is significantly larger than the kinetic 
energy term (bandwidth W^2 eV), this small difference 
in connection with the smaller bandwidth of the ai^ or- 
bitals will favor a localization of the electrons in the ai^ 
orbital for energetic reasons. 



C. Paramagnetic susceptibility: Competition of 
two exchange processes 

The LDA+DMFT(QMC) result that roughly one elec- 
tron is localized in the ai^ orbital has several immedi- 
ate consequences, which can be tested with experimental 
findings. A rather direct consequence is that the electron 
localized on the ai^ orbital leads to a local moment cor- 
responding to S = 1/2 per Vanadium atom and thus a 
Curie-like susceptibility, consistent with experiment. The 
remaining 0.5 electrons per V ion in the metallic e^ band 
will lead to a small and temperature independent Pauli 
contribution to the susceptibility. 

From an experimental point of view, the mag- 
netic properties of LiV204 pose several puz- 

x{T) = C / {T — 6) -\- xo in the temperature range 50- 
1000 K. The best fit to the experimental data is obtained 
under the assumption that the magnetic susceptibility 
x(T) is the sum of a Curie- Weiss and a T-independent 
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part xo, which contains Pauh paramagnetic, core 
diamagnetic and orbit^ Van Vleck contributions to 
the total susceptibihtyJl3 If ah V-Sd electrons (1.5 per 
V ion) would equally contribute to the formation of a 
local moment, one would expect a mixture of S = 1 and 
S = 1/2 vanadium ions and consequently an effective 
paramagnetic moment 2.34. If only one localized elec- 
tron in the ai^ orbital contributes to the Curie constant 
the effective paramagnetic moment is 1.73. However, the 
experimentally observed values of the Curie constants, 
as discussed above, in general are close to the latter 
value of the effective paramagnetic moment or slightly 
larger. Under the assumption that the local spin S 
is exactly 1/2 one gets a slightly enchanced g-factor 
from experiment. Whereas the negative Curie- Weiss 
temperature in the range ^ = —20 ... — 60 K indicates 
antiferromagnetic interactions between local S = 1/2 
moments, the large ^- value points towards ferromagnetic 
interactiojJiS|— between local moments and conduction 
electrons.BEl 

A similar ambiguity is revealed by the results of ^Li 
NMR measurements. The spin-lattice relaxation rate 
l/Ti(T) shows a broad maximum around 30-50 K and 
becomes almost temperature independent above 400 K. 
The high temperature data (T > 50 K) have been in- 
terpreted rather as an indication of localized magnetic 
moments,ll2l or, on the contrary, as characteristics of 
an itinerant, electron system close to a ferromagnetic 
instability.t^ 

The debate on the relevant magnetic interactions 
was continued after the spin fluctuation spectrum of 
LiV204 was determined by means of quasielastic neu- 
tron scattering .113 These neutron data showed a tran- 
sition from ferromagnetic correlations at elevated tem- 
peratures (T > 40 K) to antiferromagnetic spin fluctua- 
tions with a wave vector q ^ 0.7 at low temperatures 
(T < 40 K), as shown in Fig. ^. However, a subsequent 
inelastic neutron scattering study on LiV2 04 merely re- 
vealed the continuous evolution of antiferromagnetic fluc- 
tuations out of a high-temperature paramagnetic state 
without indications rof ferromagnetic interactions at ele- 
vated temperatures.E3 

In order to clarify these inconsistencies, most re- 
cently, a polarized neutron scattering study pOn a new 
set of LiV2 04 samples has been performed.o Full 3- 
directional polarization analysis allows for an unambigu- 
ous separation of the nuclear, magnetic and spin inco- 
herent cross-section, respectively. The measured mag- 
netic cross-section has been fitted and the corresponding 
real-space spin correlations extracted by employing the 
reverse Monte Carlo method. The data show a temper- 
ature induced cross-over from purely ferromagnetic next 
nearest neighbor spin correlations at high temperatures 
(T > 40 K) to a coexistence of ferromagnetic (nearest 
neighbor) and antiferromagnetic (second nearest neigh- 
bor) correlations at low temperatures. The correspond- 
ing oscillatory behavior of the real-space spin correla- 
tions at low temperature can be described purely phe- 




Q (A") 



FIG. 9: Q-dependence of the energy-integrated quasielastic 
scattering intensity of LiV2 04 at 4, 15, and 80 K. For com- 
parison, the inset shows the Q-dependence according to the 
square of the neutron magnetic form factor of the two Vana- 
dium configurations possible in LiV2 04. At elevated temper- 
atures, the much stronger reduction of the scattering intensity 
upon increasing Q points towards ferromagnetic spin correla- 
tions. The maximum around Q — 0.7 at 4 K indicates 
antiferromagnetic fluctuations. Dashe^nor dotted lines are 
guides to the eye only (taken from Ref.E3). 



nomenologically assuming a RKKY-type interaction with 
S(S+il=0.2 and kp = 0.45 A"^ for a simple parabolic 
bando Let us estimate the band filling per Vanadium 
ion that this particular value for kp corresponds to. For 
a parabolic band we have 



1 



k<kF ,cr 



a y 
2^) 



where a = 8.227A is the lattice constant. For kp = 
0.45A-1 we obtain n ^ 1.714 per cubic unit cell, i.e., a 
conduction band filling of 0.43 per Vanadium ion, since 
there are four Vanadium ions per unit cell. Note that 
this value perfectly coincides with the filling of the e^ 
band as obtained within our LDA/DMFT calculation, 
thus giving further experimental confirmation of our pro- 
posed picture of a separation of the Vanadium d-states 
into a strongly correlated ai^ orbital with one localized 
electron and weakly localized e^ orbitals which form a 
metallic band filled with ^ 0.5 electrons. 

The evolution of antiferromagnetic spin correlations 
with q ^ 0.7 A~^ at low temperatures is now experimen- 
tally well established. However, the temperature-induced 
cross-over from ferromagnetic to RKKY-like spin corre- 
lations is still unclear at present. 

We performed LDA+DMFT(QMC) calculations of the 
high-temperature magnetic susceptibility for a set of dif- 
ferent inverse temperatures /3=5, 7, 9, 10, 11, 12, 13 eV~^ 
and magnetic fields 5=0.005, 0.01, 0.02 eV. The results 
are shown as triangles for a Hund's exchange coupling 
J=0.8 eV and J=0 in Fig. 10 together with a fit to a 
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FIG. 10: Inverse Curie- Weiss spin susceptibility x ^(T) 
for LiV204 obtained by LDA+DMFT(QMC) calculation at 
[7=3.0 eV, J=0 (diamonds) an^J=0.8 eV (triangles), in com- 
parison to experimental datatZI (circles). For i) J=0 the 
Curie- Weiss constant is Q ?^-450 K and the effective para- 
magnetic moment is pe//=1.30; for ii) J=0.8 eV we obtain 
Q ?^640 K and pe//=1.65. Solid lines are least squares fits to 
the LDA+DMFT(QMC) data. 



Curie- Weiss la^ 



X{T) = Na 



2 2 
3kB{T-0) 



The resulting value of the Curie- Weiss temperature is 
(9 ^ 640 K for J = 0.8 eV, the fitted Curie- Weiss 
constant C gives an effective magnetic moment p^jj = 
g'^S{S + 1) = 3.28, which corresponds to an effective 
paramagnetic moment Peff = 1.65. The deviation from 
Peff = 1.73 for a spin 5* = 1/2 can be accounted for if 
one recalls that the occupancy of the ai^ orbital is not 
exactly 1 but 0.88. Such a reduced occupancy should 
ideally lead to a decrease of Pe// from 1.73 to 1.56. The 
slightly larger (calculated) value Peff = 1.65 is due to 
contributions from the eg states. For the calculation with 
J=0, where ai^ and e^ orbitals do not couple via a local 
exchange interaction the value of the effective paramag- 
netic moment is 1.3, which agrees with the expectations 
if one takes into account that for J=0 the occupancy of 
the diig orbital is 0.75, and hence Pe// = 1.73x0.75 = 1.3. 

While our value of Peff is in good agreement with 
known experimental data, the large ferromagnetic Curie- 
Weiss temperature of about 640K is in contrast with 
experiment— I A similar result was obtained by Anisi- 
mov et al.]^ where a rather strong effective ferromag- 
netic intersite exchange parameter J^ex = 530 K (which 
is the sum of direct and double exchanges) was obtained. 
A ferromagnetic exchange coupling between local mo- 
ments in LiV204 can be readily understood in the dou- 
ble exchange picture. The presence of two types of d- 
electrons, localized ones forming local moments and delo- 
calized ones producing a partially filled, relatively broad 



band, is a necessary ingredient for the double exchange 
mechanism, resulting in a strong ferromagnetic coupling 
between local moments. This requirement is fulfilled in 
the case of LiV2 04 with one electron localized in the ai^ 
orbital and a 1/8 filled broad e^ band. 

A second important condition for the double exchange 
mechanism is Hund's intra-atomic exchange, i.e., the 
ferromagnetic exchange interaction between electrons 
within the t2g states. If one switches off this intra- 
atomic exchange interaction, the double exchange mech- 
anism will be switched off, too. This becomes appar- 
ent from the result of our calculations for x~^(T) with 
Hund's exchange J equal to zero (see Fig. |l^, diamonds). 
Again, we obtain a Curie- Weiss like behavior with an 
Peff = 1.3. This time, however, the Curie- Weiss tempera- 
ture Odir ^ —450 K is negative^ i.e., pointing to an ef- 
fective antiferromagnetic exchange. In the absence of 
intra-atomic exchange and ai^-e^ hybridization, the only 
contribution to an exchange coupling can arise from the 
direct hybridization between the ai^ electrons. Obvi- 
ously, this will result in the observed antiferromagnetic 
exchange coupling. 

These results clearly show that, neglecting ai^-e^ hy- 
bridization, one has to expect a subtle competition in 
LiV204 between antiferromagnetic direct exchange re- 
sulting from ai^-ai^ hybridization and ferromagnetic 
double exchange from e^-e^ hybridization. Obviously, 
for the present set of Coulomb parameters the ferromag- 
netic contribution wins in our DMFT calculations. Note, 
that this result is a direct consequence of the value of J 
used in our calculations. There are also available experi- 
mental results from high-energy spectroscopgqwhich lead 
one to expect a value of J = 0.65 ... 0.7 eV O instead of 
J = 0.8 eV used in our calculation. We found that such a 
smaller value of J do not change our results significantly. 

There is, however, an important part missing in our 
present calculation. Since we use the LDA DOS and not 
the full Hamiltonian in the DMFT caclulations, the aig- 
Cg hybridization is neglected completely. As has been 
pointed out by Anisimov et a/.H, this hybridization can 
give rise to an antiferromagnetic Kondo coupling between 
aig and orbitals. The value of this additional exchange 
coupling was estimated to be Jx=-630 Ko Thus, in 
LiV204 there are possibly three important exchange in- 
teractions present between the ai^ and e^ electrons of the 
V3d shell: ferromagnetic double exchange (~ 1090 K), 
antiferromagnetic direct exchange (?^-450 K) and antifer- 
romagnetic Kondo exchange (~-630 K). All these inter- 
actions effectively cancel, eventually leading to a small 
Curie- Weiss temperature of the order of K in accor- 
dance with experiment. 

Finally, one puzzle remains, namely that the data from 
neutron scattering experiments shpw-ierromagnetic spin 
fluctuations at high temper at ures.t§LJ 

The competing exchange terms described above may 
also explain the change from antiferromagnetic to ferro- 
magnetic spin fluctuations at T 40 K found in neutron 
scattering experiments. This is about the same energy 



10 



scale as the Kondo temperature. Thus, one might argue 
that, while below Tk the combined antiferromagnetic di- 
rect and Kondo exchanges are stronger than the ferro- 
magnetic double exchange, this changes above Tk where 
the Kondo effect becomes ineffective, such that the fer- 
romagnetic exchange prevails. 

IV. CONCLUSION 

We investigated the effect of Coulomb correlations on 
the electronic structure, orbital state and magnetic prop- 
erties of LiV204. The analysis of the non-interacting 
partial densities of state obtained by standard LDA cal- 
culations shows that, while the trigonal splitting of the 
t2g states into ai^ and e^ orbitals is not strong enough 
to produce separate bands, it leads to a significant dif- 
ference in the effective bandwidths and average energy of 
the trigonal orbitals. The LDA+DMFT(QMC) calcula- 
tions gave orbital occupancies and spectra which indicate 
a nearly complete localization of one electron out of 1.5 
3(i electrons per V ion in the ai^ orbital, while the e^ 
orbitals form a relatively broad partially filled metallic 
band. The calculated temperature dependence of the 
paramagnetic susceptibility corresponds to the experi- 
mentally observed Curie- Weiss law and gives an effec- 



tive paramagnetic moment j9e// = l-65 in agreement with 
experimental data. The experimentally observed small 
value of the Curie- Weiss temperature, formerly a puzzle, 
is supposed to be the result of a competition between 
three different contributions to the effective exchange in- 
teraction between ai^ and e^ electrons in the V3d shell: 
ferromagnetic double exchange, antiferromagnetic direct 
exchange and antiferromagnetic Kondo exchange. The 
present calculations show a dominance of the first, but 
including the hybridization between aig and e^, an al- 
most cancellation of all three becomes possible leading 
to the experimentally observed small residual antiferro- 
magnetic Curie- Weiss temperature. 



Acknowledgments 

We thank A. Oles, A. Loidl, D. Singh, M. Korotin, 
P. Horsch and V. Eyert for very useful discussions 
and A. Sandvik for making available his maximum en- 
tropy program. This work was supported in part by 
the Deutsche Forschungsgemeinschaft through Sonder- 
forschungsbereich 484 (DV, TP, GK, AK, IN) and the 
Russian Foundation for Basic Research through grants 
RFFI-01-02-17063 (VA,IN) and RFFI-02-02-06162 (IN). 



^ K. Andres, J.E. Graebner, and H.R. Ott, Phys. Rev. Lett. 

35, 1779 (1975). 
2 G.R. Stewart, Rev. Mod. Phys. 56, 755 (1984). 
^ G. Aeppli and Z. Fisk, Comments Condens. Matter Phys. 

16, 155 (1992). 

^ S. Kondo, D.C. Johnston, C.A. Swenson, F. Borsa, 
A.V. Mahajan, L.L. Miller, T. Gu, A.I. Goldman, 
M.B. Maple, D.A. Gajewski, E.J. Freeman, N.R. Dilley, 
R.P. Dickey, J. Merrin, K. Kojima, G.M. Luke, Y.J. Ue- 
mura, O. Chmaissen, and J.D. Jorgensen, Phys. Rev. Lett. 
78, 3729 (1997). 

^ R. Ballon, E. Lelievre-Berna, and B. Fak Phys. Rev. Lett. 
76, 2125 (1996), and references therein. 

^ S.A. Carter, T.F. Rosenbaum, P. Metcalf, J.M. Honig, and 
J. Spalek, Phys. Rev. B 48, 16841 (1993). 
C. Urano, M. Nohara, S. Kondo, F. Sakai, H. Takagi, 
T. Shiraki, and T. Okubo, Phys. Rev. Lett. 85, 1052 
(2000). 

^ K. Kadowaki and S. B. Woods, Sol. St. Com. 58, 507 
(1986). 

^ D.C. Johnston, C.A. Swenson, and S. Kondo, Phys. Rev. 
B 59, 2627 (1999). 

T. Hayakawa, D. Shimada and N. Tsuda, J. Phys. Soc. 
Jpn. 58, 2867 (1989). 

S. Kondo, D.C. Johnston, and L.L. Miller, Phys. Rev. B 
59, 2609 (1999). 

M. Onoda, H. Imai, Y. Amako, and H. Nagasawa, Phys. 
Rev. B 56, 3760 (1997). 

M. Lohmann, J. Hemberger, M. Nicklas, H.A. Krug von 
Nidda, A. Loidl, M. Klemm, G. Obermeier, and S. Horn, 
Physica B 259-261, 963 (1999). 



A. Fujimori, K. Kawakami and N. Tsuda, Phys. Rev. B 
38, 7889 (1988). 

N. Fujiwara, H. Yasuoka, and Y. Ueda, Phys. Rev. B 57, 
3539 (1998). 

N. Fujiwara, H. Yasuoka, and Y. Ueda, Phys. Rev. B 59, 
6294 (1999). 

A.V. Mahajan, R. Sala, E. Lee, F. Borsa, S. Kondo, and 

D.C. Johnston, Phys. Rev. B 57, 8890 (2000). 

A. Krimmel, A. Loidl, M. Klemm, H. Schober, and S. Horn, 

Phys. Rev. Lett. 82, 2919 (1999); Phys. Rev. B 61, 12578 

(1998); Physica B 276-278 766 (2000). 

W. Trinkl and A. Loidl, M. Klemm, and S. Horn, Phys. 

Rev. B 62, 8915 (2000). 

S.-H. Lee, Y. Qiu, C. Broholm, Y. Ueda, and J.J. Rush, 
Phys. Rev. Lett. 86, 5554 (2001) 

K. Fujiwara, H. Yoshioka, K. Miyoshi, J. Takeuchi, 
T.C. Kobayashi, and K. Amaya, Physica B 312-313, 913 
(2002). 

D.C. Johnston, |cond-mat/9910404 . 

F. Takagi, K. Kawakami, L Maekawa, Y. Sakai, and 

N. Tsuda, J. Phys. Soc. Jpn. 56, 444 (1987). 

Y. Nakajima, Y. Amamiya, K. Ohnishi, L Terasaki, 

A. Maeda, and K. Uchinokura, Physica C 185-189, 719 
(1991). 

Y. Ueda, N. Fujiwara, and H. Yasuoka, J. Phys. Soc. Jpn. 
66, 778 (1997). 

H. Kessler and M.J. Sienko, J. Chem. Phys. 55, 5414 
(1971). 

B. L. Chamberland and T.A. Hewston, Solid State Comm. 
Vol.58, No. 10, pp.693-695 (1986). 

M. Brando, N. Biittgen, V. Fritsch, J. Hemberger, H. Kaps, 



11 



H.-A. Kmg von Nidda,M. Nicklas, K. Pucher, W. Trinkl, 
A. Loidl, E.W. Scheldt, M. Klemm, and S. Horn, Eur. 
Phys. J. B25, 289-298 (2002). 

K. Miyoshi, M. Ihara, K. Fujiwara, and J. Takeuchi, Phys. 
Rev. B 65, 092414 (2002). 
^° D. C. Johnston, T. Ami, F. Borsa, M. K. Crawford, J. A. 
Fernandez-Baca, K. H. Kim, R. L. Harlow, A. V. Mahajan, 
L. L. Miller, M. A. Subramanian, D. R. Torgeson and Z. 
R. Wang, Springer Ser. Solid State Sci. 119, 241 (1995) 
W. Kohn, L.J. Sham, Phys. Rev. A - Gen. Phys. 140, 
1133 (1965); L.J. Sham, W. Kohn, Phys. Rev. 145 N 2, 
561 (1966). 

L. Hedin and B. Lunqvist, J. Pys. C: Solid State Phys. 4, 
2064 (1971); U. von Barth and L. Hedin, J. Pys. C: Solid 
State Phys. 5, 1629 (1972). 

V.L Anisimov, M.A. Korotin, M. Zolfl, T. Pruschke, 
K. Le Hur, and T.M. Rice, Phys. Rev. Lett. 83, 364 (1999). 
O.K. Andersen, Phys. Rev. B 12, 3060 (1975); O. Gun- 
narsson, O. Jepsen, and O.K. Andersen, Phys. Rev. B 27, 
7144 (1983). 

V.L Anisimov, J. Zaanen, and O.K. Andersen, Phys. Rev. 
B 44, 943 (1991);J. Phys.: Condens Matter 9, 767 (1997). 
2^ V. Eyert, K.-H. Hock, S. Horn, A. Loidl, and P.S. Risebor- 
ough, Europhys. Lett. 46, 762 (1999). 
A.R. Williams, J. Kiibler, and CD. Jr. Gellat, Phys. Rev. 
B 19, 6094 (1979). 

J. Matsuno, A. Fujimori, and L.F. Mattheiss, Phys. Rev. B 
60, 1607 (1999). 

L.F. Mattheiss and D.R. Hamann, Phys. Rev. B 33, 823 
(1986). 

O.K. Andersen, Phys. Rev. B 12, 3060 (1975). 

L.F. Mattheiss and W. Weber, Phys. Rev. B 25, 2258 

(1982). 

D.J. Singh, P. Blaha K. Schwarz, and LL Mazin, Phys.Rev. 
B 60, 16359 (1999). 

D.J. Singh, Planewaves, Pseudopotentials and the LAPW 
Method (Kluwer Academic, Boston, 1994). 
P. Blaha, K. Schwartz, and J. Luitz, Vienna University 
of Technology, 1997 [Improved and updated Unix version 
of the WIEN code, published by P. Blaha, K. Schwartz, 
P. Sorantin, and S.B. Trickey, Comput. Phys. Commun. 
59, 399 (1990)]. 

S.H. Wei and H. Krakauer, Phys. Rev. Lett. 55, 1200 

(1985); D. Singh, Phys. Rev. B 43, 6388 (1991). 

CM. Varma, Phys. Rev. B 60, R6973 (1999). 

H. Kusunose, S. Yotsuhashi, and K. Miyake, Phys. Rev. B 

62, 4403 (2000). 

J. Hopkinson and P. Coleman, |cond-mat / 0203288| . 
J. Hopkinson and P. Coleman, Physica B 312, 711 (2002). 
S. Fujimoto, Phys. Rev. B 65, 155108 (2002). 
S. Burdin, D.R. Grempel, and A. Georges, |cond- 



52 
53 



at/0107284 



B. Renter and J. Jaskowsky, Angew. Chem. 72, 209 (1960). 
H. Wada, M. Shiga, and Y. Nakamura, Physica B 161, 197 
(1989). 

O. Chmaissem and J.D. Jorgensen, S. Kondo, and 
D.C. Johnston, Phys. Rev. Lett. 79, 4866 (1997). 
V. Eyert and K.-H. Hock, Phys. Rev. B 57, 12727 (1998). 
K. Terakura, T. Oguchi, A.R. Williams, and J. Kiibler, 
Phys. Rev. B 30, 4734 (1984). 

W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989). 

D. Vollhardt, in Correlated Electron Systems, edited by 



V. J. Emery, World Scientific, Singapore, 1993, p. 57. 
Th. Pruschke, M. Jarrell, and J. K. Freericks, Adv. in Phys. 
44, 187 (1995). 

A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

V. L Anisimov, A. L Poteryaev, M. A. Korotin, A. O. 
Anokhin, and G. Kotliar, J. Phys. Cond. Matter 9, 7359 
(1997). 

A. Lichtenstein and M. L Katsnelson, Phys. Rev. B 57, 
6884 (1998). 

K. Held, LA. Nekrasov, N. Bliimer, V.L Anisimov, and 
D. Vollhardt, Int. J. Mod. Phys. B 15, 2611 (2001); K. 
Held, LA. Nekrasov, G. Keller, V. Eyert, N. Bliimer, A.K. 
McMahan, R.T. Scalettar. T. Prusc hke. V.L Anisimov, 
and D. Vollhardt, |cond-mat/0112079| (Pubhshed in Quan- 
tum Simulations of Complex Many-Body Systems: From 
Theory to Algorithms, eds. J. Grotendorst, D. Marks, and 
A. Muramatsu, NIC Series Volume 10 (NIC Directors, 
Forschunszentrum Jiilich, 2002) p. 175-209. (ISBN 3-00- 
009057-6) 

J. E. Hirsch and R. M. Eye, Phys. Rev. Lett. 56, 2521 
(1986); M. Jarrell, Phys. Rev. Lett. 69, 168 (1992); M. 
Rozenberg, X. Y. Zhang, and G. Kotliar, Phys. Rev. Lett. 
69, 1236 (1992); A. Georges and W. Krauth, Phys. Rev. 
Lett. 69, 1240 (1992); M. Jarrell in Numerical Methods 
for lattice Quantum Many- Body Problems, edited by D. 
Scalapino, rAddison Wesley, 1997. For review of QMC for 
DMFT seeH 

For a review on the maximum entropy method see M. Jar- 
rell and J. E. Gubernatis, Physics Reports 269, 133 (1996). 
I. A. Nekrasov, K. Held, N. Bliimer, A. I. Poteryaev, V. 
I. Anisimov, and D. Vollhardt, Euro. Phys. J. B 18, 55 
(2000). 

K. Held, G. Keller, V. Eyert, V. I. Anisimov, and D. Voll- 
hardt, Phys. Rev. Lett. 86, 5345 (2001). 
M. B. Zolfl, Th. Pruschke, J. Keller, A. I. Poteryaev, I. 
A. Nekrasov, and V. I. Anisimov, Phys. Rev. B 61, 12810 
(2000). 

A. Liebsch and A. Lichtenstein, Phys. Rev. Lett. 84, 1591 
(2000). 

Ph. Lambin and J. P. Vigneron, Phys. Rev. B 29, 3430 
(1984). 

O. Gunnarsson, O.K. Andersen, O. Jepsen, and J. Zaanen, 
Phys. Rev. B 39, 1708 (1989); V.L Anisimov and O. Gun- 
narsson, Phys. Rev. B 43, 7570 (1991). 
Recently the Curie- Weiss law was derived analitically 
in dynamical mean-field theory for the Hubbard model 
(K. Byczuk and D. Vollhardt, Phys. Rev. B 65, 134433 
(2002)). 

J. Zaanen and G. A. Sawatzky extracted the value of 
J=0.64 eV obtained from Racah parameters for V ions 
(J. Solid State Chem. 88, 8 (1990)) using the procedure, 
described in J. S. Griffith, The Theory of Transition Metal 
Ions (Cambridge University Press, Cambridge, 1971). T. 
Mizokawa and A. Fujimori deduced from the cluster- model 
analyses of the photoemission spectra J=0.68 eV (Phys. 
Rev. B 54, 5368-5380 (1996)). 

A. Krimmel, A. P. Murani, J. R. Stewart, A. Ibarra-Palos, 
P. Strobel and A. Loidl, Acta Physica Polonica, in print 
(Proceedings of the International Conference on Strongly 
Correlated Electron Systems, Krakov 2002) 



